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ABSTRACT 


We provide an algorithm for enumerating near-minimum weight s-t cuts in 
directed and undirected graphs, with applications to network interdiction and network 
reliability. “Near-minimum” means within a factor of 1+8 of the minimum for some 
8 > 0. The algorithm is based on recursive inclusion and exclusion of edges in locally 
minimum-weight cuts identified with a maximum flow algorithm. We prove a 
polynomial-time complexity result when 8 = 0, and for 8 > 0 we demonstrate good 
empirical efficiency. The algorithm is programmed in Java, run on a 733 MHz Pentium 
III computer with 128 megab 3 ^es of memory, and tested on a number of graphs. For 
example, all 274,550 near-minimum cuts within 10% of the minimum weight can be 
obtained in 74 seconds for a 627 vertex 2,450 edge unweighted graph. All 20,806 near¬ 
minimum cuts within 20% of minimum can be enumerated in 61 seconds on the same 
graph with weights being uniformly distributed integers in the range [1,10]. 
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EXECUTIVE SUMMARY 


We consider the problem of enumerating all near-minimum weight s-t cuts in a 
directed graph G=iV,E). A minimal s-t cut is a minimal set of edges whose removal 
results in disconnection of the “source vertex” s and “sink vertex” t. A cut is “near 
minimum” if its weight is within a multiplicative factor of l+£ times the minimum cut 
weight for some £ > 0. 

This thesis is motivated by network interdiction problems in which a decision 
maker selects edges to interdict based on a primary criterion, e.g., resource consumption, 
and some secondary criteria, e.g., collateral damage, exposure to enemy fire, etc. 
However, the results may also be useful for evaluating graph reliability and other 
problems. 

Our enumeration algorithm is based on a recursive “inclusion-exclusion” method. 
The algorithm identifies a (minimal) minimum-weight s-t cut by exploiting the duality 
between max flow and minimum cut, and then systematically partitions the space of 
minimal (and possibly some non-minimal) cuts by attempting to include and exclude the 
edges in that cut. The quasi-exclusion of an edge is accomplished by setting its weight to 
infinity, so that the edge will not be identified as a part of any finite-weight cut in 
subsequent solutions. We quasi-include an edge (m,v) by implicitly adding the infinite 
weight edges {s,u) and (v,0 to the graph so that every finite-weight cut in G must now 
contain (m,v). The weight of the locally minimum cut is monotonically non-decreasing, 
so the recursive algorithm can backtrack as soon as the weight of a locally minimum cut 
exceeds (l+£)W)v//Ar, where Wmin is the minimum weight of a cut in G. 

xiii 


The quasi-inclusion technique can, unfortunately, lead to enumeration of non- 
minimal cuts, i.e., cuts that contain more edges than necessary to disconnect s from t. We 
identify non-minimal cuts by performing two searches beginning from s and t and not 
traversing cut edges: If any endpoint of a cut edge is not reached, then the cut is non- 
minimal. The algorithm simply ignores non-minimal cuts and continues. 

The algorithm is correct because the sequence of weights of (locally) minimum 
cuts is monotonic non-decreasing as we descend along any path in the inclusion- 
exclusion enumeration tree and because we never eliminate or duplicate any minimi, 
near-minimum cut during enumeration. 

The complexity of the algorithm for finding only minimum cuts (when 8 = 0) is 
0(/(lVI,l£:i)-l-IVIIi?IICo(G)i) where /(iyi,l£:i) is the complexity of solving a 

maximum flow problem on G = (y,E) from scratch and Co(G) is the set of minimum cuts 
in G. Unfortunately, this polynomial complexity result is not valid for near-minimum cut 
enumeration, because we cannot properly bound the number of non-minimal cuts 
identified during the enumeration. Thus the worst-case complexity of the algorithm for 
near-minimum cut enumeration remains unknown when e > 0. 

The enumeration algorithm is implemented in Java (JDK 1.2.2) and tested using a 
733 MHz Pentium HI computer with 128 megabytes of RAM operating under Microsoft 
Windows 98 SE. For example, in a 25 by 25 grid graph with 627 vertices, 2,450 edges 
and unit weights, all 24 minimum-weight cuts (8 = 0) are enumerated in 0.22 seconds, all 
1,128 cuts with 8 = 0.05 are enumerated in 3.91 seconds, and all 3,621,978 cuts with 
8 = 0.15 are enumerated in 973.29 seconds. 
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Excerpt from Walden 

I went to the woods because I wished to live deliberately^ to front only the essential facts of life, and see if 1 
could not learn what it had to teach, and not, when I came to die, discover that / had not lived. I did not wish 
to live what was not life, living is so dear, nor did I wish to practice resignation, unless it was quite necessary. 
I wanted to live deep and suck all the marrow of life, to live so sturdily and Spartan-like as to put to rout all 
that was not life, to cut a broad swath and shave close, to drive life into a corner, and reduce it to its lowest 
terms, and if it proved to be mean, why then to get the whole and genuine meanness of it, and publish its 
meanness to the world; or if it were sublime, to know it by experience, and be able to give a true account of it 
in my next excursion. For most men, it appears to me, are in a strange uncertainty about it, whether it is of the 
devil or of God, and have somewhat hastily concluded that it is the chief end of man here to glorify God and 
enjoy him forever ." 

- Henry David Thoreau 
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I. 


INTRODUCTION 


Enumeration of various types of minimal cuts in graphs has been well studied, 
and both the operations research and reliability engineering communities have presented 
efficient enumeration algorithms. In this thesis, we deal with a particular type of cut that 
has not received the same attention as other types, specifically, near-minimum-weight 
minimal s-t cuts. 

We develop and implement an empirically efficient algorithm to enumerate all 
near-minimum-weight minimal s-t cuts in a directed or undirected graph. (Note: “cuts” 
will be used synonymously with “minimal cuts” from here on unless stated otherwise. 
Also, note that “minimal” refers to the set-theoretic property that a set of edges that 
disconnects all s-t paths is minimal, and “minimum” refers to the total weight of a cut.) 

A. NEAR-MINIMUM CUTS AND THEIR APPLICATIONS 

The network to be considered, primarily, is a directed graph G = {V, E), with 
positive edge weights and two special, distinct vertices, a source s and a sink t. An s-t cut 
is a minimal set of edges whose removal breaks all directed paths from s to t. We define 
the problem of finding an s-t cut of minimum weight among all possible s-t cuts in G as 
the minimum s-t cut problem (MCP-st). We also define two extensions of MCP-st, the 
problem of enumerating all minimum s-t cuts in G (AMCP-st) and the problem of 
enumerating all near-minimum weight s-t cuts (ANMCP-st) whose weight is within a 
factor of 1-i-e of the minimum s-t cut weight in the given network. The last problem is the 
focus of this thesis. 
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The analogs of AMCP-st and ANMCP-st in undirected weighted graphs G are 
easily solved using the techniques we develop, also. An s-t cut in an undirected network 
is simply a minimal set of edges C whose deletion breaks all (undirected) s-t paths. If we 
replace each undirected edge by two directed, anti-parallel edges, each with the weight of 
the original undirected edge—this is a standard transformation, of course—each s-t cut in 
the new directed graph corresponds in a one-to-one fashion with an s-t cut in the original 
graph. Thus, an efficient technique to enumerate s-t cuts in a directed graphs leads to an 
efficient technique to enumerate s-t cuts in an undirected graph. 

Another type of cut can be defined in an undirected graph G which we call a 
general cut, this is simply any minimal disconnecting set of edges in G. The problems of 
finding and enumerating general cuts are somewhat related to our problems and will be 
discussed briefly, so we make these definitions: (a) The problem of finding a minimum- 
weight general cut is the general min-cut problem (MCP), (b) the problem of 
enumerating all minimum-weight general cuts is the general all min-cuts problem 
(AMCP), and (c) the problem of enumerating all near-minimum-weight general cuts is 
the general all near-min-cuts problem (ANMCP). 

In the remainder of this thesis, we will use the following notation and define more 
as needed. For s-t cuts, we denote a minimum cut as Co and a near-minimum cut as 
with their respective cardinalities ICol and IC^ I. We also let Co(G) and C^(G) denote the 

set of mi nimum and near-minimum cuts in G, respectively. 

One application of ANMCP-st arises in network interdiction problems (Wood 

1993). A network user (adversary) attempts to “communicate” between source s and 

s ink r in a directed network while the interdictor (decision maker), using limited 
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resources (aerial sorties, cruise missiles, etc.), tries to break all paths between 5 and t. By 
treating the amount of resource required to destroy an edge as its weight, a decision 
maker can use a maximum flow minimum-cut algorithm to identify a minimum-weight 
cut, i.e., minimum-resource cut, to disrupt the communication between s and t. 
(Hereafter, following common practice, the term max is substituted for maximum and 
min for minimum.) 

There may be secondary criteria that the decision maker wishes to consider when 
determining the best interdiction plan, e.g., collateral damage, risk to attacking forces, 
etc. In this case, near-optimal solutions with respect to the primary criterion can be 
obtained by solving ANMCP-st; then those solutions can be evaluated against the 
secondary criteria for suitability. One of those near-optimal “good solutions” might 
produce more desirable results than an “optimal solution” obtained by solving MCP-st or 
AMCP-st (Boyle 1998, Gibbons 2000.) 

Another application of ANMCP-st arises in assessing the reliability and 
connectivity of networks (Colboum 1987, Provan and Ball 1983). Consider a 
probabilistic graph G with two specified vertices s and /, in which each edge fails 
independently with probability p. The s-t connectedness (or two-terminal reliability in 
the undirected case) of the graph is the probability that there exists a working s-t path 

l£1 

in G. This probability can be computed as 1- ^ p* (1 - p)'^'“*, where kmin is the 

cardinality of the minimum cut and \ is the number of edge sets of cardinality k whose 
deletion disconnects s from t. When p is sufficiently small, this value can be 
approximated by l-y^jj^j^p*"'"(l-p)'^'”*'"“, which equals l-ICo(G)l p'‘^°'(l-when 
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edge weights are all defined to be 1 (i.e., weight of a cut corresponds to the number of 
edges in the cut) (Colboum 1987, pp. 53-54). This approximation (which is also an upper 
bound) can be made more accurate by including terms from the original formula for 
k^^<k<k', where k' is “not very large.” (Intuitively, only cuts with small cardinality 
can have a significant probability of failing, and thus disconnecting 5 and t; Karger 1999) 
When k = itmin, \ counts only minimal cuts, but when k > k^in, non-minimal cuts must 
be counted in addition to minimal ones. However, non-minimal cuts of cardinality k are 
easily counted given the minimal cuts of cardinality ^-1, k-2, etc. Thus, enumeration of 
minimum and near-minimum cuts is essential for accurately estimating s-t connectedness 
ofG. 

B. BACKGROUND 

As a result of their widespread applicability in network reliability and in many 
other combinatorial problems (Picard and Queyranne 1982), AMCP, AMCP-st and the 
problem of enumerating all minimal s-t cuts have been intensively studied, but ANMCP- 
st has not received the same attention. 

One bmte-force approach for ANMCP-st would be to find all minimal s-t cuts and 
then eliminate the ones that are not minimum or near-minimum. All minimal s-t cuts of a 
graph can be enumerated in time that is polynomial in the size of G and in the number of 
cuts enumerated (Tsukiyama et al., 1980, Abel and Bicker 1982, Karzanov and Timofeev 
1986, Shier and Whited 1986, Sung and Yoo 1992, Ahmad 1990, Prasad et al., 1992, 
Nahman 1995, Patvardhan et al., 1995 and Fard and Lee 1999). Unfortunately, this 
cannot lead to an efficient general approach for AMCP-st or ANMCP-st because the 
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number of cuts in a graph may be exponential in the size of that graph while the number 
of minimum and near-minimum cuts may be polynomial. For instance, if G is a complete 
directed graph with edge weights of 1, the total number of minimal cuts is the 

number of minimum cuts is 2 and the number of (minimal) cuts of the next largest size is 
2(IVI-2). 

All minimum s-t cuts (AMCP-st) can be enumerated efficiently. Picard and 
Que 5 a:anne (1980) present a characterization of all minimum s-t cuts to solve AMCP-st 
by exploiting “max flow-min cut” duality. They find a maximum s-t flow in G, define a 
binary relation (Picard and Queyranne 1980, Theorem 1) on the vertices of the resulting 
“residual graph” and show how to find all minimum s-t cuts by efficiently generating all 
closures of this graph. Gusfield and Naor (1993) and Curet (1999) give efficient 
algorithms for AMCP-st using this idea. Ball and Provan (1983) identify all “ 5 -directed 
minimum cuts” by solving AMCP-st between s and i, for each i e V - { 5 }, and upon 
identification of the solution to AMCP-si, they collapse vertices s and / into a single 
vertex s to find the remaining .y-directed cuts (which break all paths between s and 
T= V- { 5 }). Provan and Shier (1996) define a general paradigm to solve AMCP-st by 
extending the Picard and Queyranne approach. Their paradigm lists s-t cuts in G(IV1) 
time per cut in directed and undirected graphs. They also show that the algorithms of 
Picard and Queyranne (1980) and Ball and Provan (1983) may not be polynomial in the 
case of directed graphs. Kanevsky (1993) uses Provan and Shier’s paradigm to find all 
minimum-size “separating vertex sets” instead of edge sets. 

Before going any further in discussing algorithms to solve ANMCP and ANMCP- 

st efficiently, we must define some notation. We let n = IVl and m = lEI from here on. 
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Ramanathan and Colboum (1987) solve ANMCP by giving an efficient algorithm to 
count “almost-minimum cardinality s-t cuts.” They bound the number of cuts enumerated 
by where it > 1 is a constant that determines the cardinality of almost- 

minimum cut together with the cardinality of the minimum cut. This algorithm is 
applicable only to undirected graphs and polynomial complexity is guaranteed only if k is 
fixed. 

Karger and Stein (1996) introduce a randomized algorithm for solving ANMCP 
by using edge contraction (identify an edge that does not cross the minimum cut and 
merge its endpoints into a single new vertex without losing the minimum cut). Their 
algorithm enumerates all general cuts whose weight is within a factor a of the minimum 
cut in 0{n^ log^ n) expected time with high probability. They also derive an upper 
bound 0(n^“)on the number of cuts. Karger (2000) later improves this upper bound to 
. But this is a Monte Carlo algorithm and gives the right answer with high 
probability, not with certainty. Nagamochi et al. (1997) give a deterministic algorithm 
similar to that of Karger and Stein (1996) for solving ANMCP. They show that all 
general cuts of weight less than k times the general minimum cut weight can be 
enumerated in Oim^n + n^'^m) time. Unfortunately, both algorithms apply only to 
undirected graphs and, furthermore, they are unlikely to be extendable to enumeration 
problems involving s-t cuts (Karger and Stein 1996). 

Vazirani and Yannakakis (1992) propose a polynomial-time algorithm for solving 
ANMCP and ANMCP-st in a directed or vmdirected graph with weighted or unweighted 
edges. In their extended abstract, they give an algorithm to find all general cuts and all 
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s-t cuts in ascending order with respect to their weight. They show that the k-th minimum 
weight s-t cuts in this order can be enumerated using maximum flow 

computations, for any fixed k. They list the s-t cuts in increasing order of the weight by 
using an idea similar to that of Picard and Queyranne (1980) (shrinking the strongly 
connected components to single vertices in the residual graph) and provide a bound of 
0(n^*"') for the number of kth minimum weight cuts in an undirected graph. The 
technique of Vazirani and Yannakakis may be useful, but it is not as easy to describe as 
ours and is based on an unproven assertion. In particular, they state: “Given a partially 
specified cut, we can find with one max flow computation a minimum weight s-t cut 
consistent with it.” We leave it to others to determine the validity of their approach. 

Boyle’s algorithm for solving constrained network-interdiction problems on 
undirected graphs (Boyle 1999) can be modified to enumerate near-minimum cuts. But 
the algorithm is only applicable to planar graphs and, because the technique is based on 
the planar dual of a graph, no generalization to non-planar graphs seems likely. (A 
generalization to directed planar graphs might be possible.) Gibbons (2000) gives an 
algorithm for solving ANMCP-st in a directed and undirected graph, but he may 
enumerate a cut more than once. Empirically, the running time and number of cuts 
enumerated grow rapidly as the size of graph and e increase. 

The current approaches for solving ANMCP-st are either computationally 

prohibitive (Gibbons 2000), not truly efficient (Ramanathan and Colboum 1987), or 

restricted to specific types of graphs, e.g., undirected (Ramanathan and Colboum 1987), 

planar (Boyle 1999). There are some efficient techniques to enumerate all near-minimum 

general cuts, i.e., to solve ANMCP (Karger 2000, Nagamochi et al. 1997), but these do 
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not seem to be extendable to ANMCP-st. Only the algorithm due to Vazirani and 
Yannakakis (1992) may be efficient for solving ANMCP-st (and ANMCP), but their 
technique is complicated and based on an unproven assertion. Therefore, we propose a 
new algorithm to solve all instances of ANMCP-st efficiently, and provide computational 
results for a simple implementation of the algorithm. 

Our algorithm for solving ANMCP-st first identifies a minimum weight s-t cut by 
solving a maximum flow problem. Given a minimum cut, the algorithm recursively 
partitions the space of possible cuts by forcing inclusion and exclusion of cut edges from 
subsequent cuts (which are possible solutions) to enumerate all possible near-minimum 
cuts in G. An edge (m,v) is excluded from a cut by setting its weight to infinity. Edge 
inclusion is accomplished by implicitly adding two infinite-weight edges to G. One edge 
connects s to u and the other connects v to r so that they create a simple path from 5 to r 
via that edge. The technique for including edges, unfortunately, can cause the algorithm 
to identify non-minimal cuts, i.e., cuts containing a minimal cut as a subset. 

C. THESIS OUTLINE 

The remainder of this thesis is organized as follows: Chapter II introduces more 
definitions and notation. In Chapter III, we explain our algorithm theoretically and prove 
its correctness. Chapter IV contains computational results from testing the algorithm on 
different types of graphs with different sizes. Finally, Chapter V summarizes our findings 
and gives conclusions. 
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II. DEFINITIONS AND NOTATION 


Let G = (V, £) be an edge-weighted directed graph with a finite set V of vertices 
and a set £■ c VxF of edges. A directed edge e = (u,v) is an incoming edge to v and an 
outgoing edge from u. An undirected graph is defined similarly, except that its edges are 
unordered pairs from V (with no orientation). We may use e or (m,v) as shorthand for 
e = (m,v). We denote the number of vertices by n = IVl and the number of edges by m = 
\E\. We distinguish two separate vertices s and tin V as the source and sink, respectively. 
Edge weights are specified by the weight function w: E ^ Z*, where is the set of non¬ 
negative integers. We denote the weight of an edge as We or w(u,v) and the vector of 
edge weights as w = , w^^w^^) . 

An s-t cut in G is a set of edges whose removal disconnects s from t in G, i.e., 
breaks all (directed) paths from s to t. An s-t cut is minimal if it does not contain a subset 
which itself is an s-t cut of the graph. We will use “s-t cut” and “cut” interchangeably for 
“minimal s-t cut” in the remainder of the thesis; a cut is denoted by C. (Note: Whenever 
we need to mention cuts other than s-t cuts, we may refer such cuts as simply general 
cuts.) The weight w{C) of an s-t cut is the sum of the edge weights in the cut, that is, 
w(C) = ^ fninimum cut (min cut) is an s-t cut whose weight is minimum among 

all s-t cuts; such a cut is denoted Cq. A near-minimum cut (near-min cut) is an s-t cut 
whose weight is at most 1 + e times the minimum cut weight for a given e > 0; such a cut 
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is denoted C^. We use Co(G) to denote the set of minimum cuts, and C^(G) to denote 
the set of near-min cuts. 

A/Z( 9 w in a directed G is a function /: £ Z, where/(e) < w{e) for each edge 
e € £, and all m e 7- {5,f}. The value of flow from 

5 to r is F = ^ ^ f(s,u). A flow F* with the maximum value among all flows is called 

maximum flow between s and t. In the maximum flow problem we wish to find a flow of 
maximum value from s to t. 

As a result of the max flow min-cut theorem and its proof (e.g., Ahuja et al 
pp. 184-185), we know that h’(Co) = F*, where Co e Co(G). Given any maximum flow F* 
together with/*=(e) for every e e £, we can identify a minimum cut Co in 0{m) time. 

In describing our algorithm, we use Wmin and Wmw to denote the weight of 
minimum-cut in G and the weight of a newly identified cut (which will be a min cut in a 
modified graph), respectively. 

A rooted tree, T, is a free tree (a connected, acyclic, undirected graph) in which 
one node (called the “root” and denoted by r) is distinguished from others. A node i in T 
with root r is said to be in level (depth) I if the length of unique path P, = ( r, vi, V 2 ,.., v^.i, 
i) from root r to node i is k. Every node along path P, (except node i) is called a (proper) 
ancestor of i, and if i is ancestor of j, then j is descendant of i. Let h and i be the last two 
nodes of P„ then h is the parent of i and i is a child of h. Nodes with the same parents are 
called siblings. A rooted tree, and enumeration tree, will be used to describe the 
recursive partitioning process that we use to solve ANMCP-st. 
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III. THEORETICAL RESULTS 


In this chapter, we introduce our algorithm (Algorithm B) for solving ANMCP-st 
and AMCP-st. We first describe the generic partitioning algorithm, Algorithm Al, from 
Gibbons (2000) as the basic approach. This is modified into Algorithm A2, which is then 
approximately implemented as Algorithm B. 

A. BASIC ALGORITHMS 

Algorithm Al (Figure 1) provides a basic framework for enumerating near-min 
s-t cuts although it may be difficult or impossible to implement efficiently. It begins by 
finding a minimum cut Co and its weight w(Co). Then, the algorithm calls the procedure 
ENUMERATE which attempts to find a new minimum cut by processing the edges of 
originally identified cut such that the edges are forced into (included in) or out of 
(excluded from) any new near-min cuts. For example, let Co = {e\, e 2 , ek} be the 
edges of the initial minimum cut. Based on this cut, (G) can be partitioned as 

C, (G) = [ C^ (C) n (e,)] u [ C, (G) n ci n (e.)] u [ C, (G) n c, n ^2 n (e^)] u 
...u [C^(G) n Cl n e2 n..n (c*)] u [C^(G) n ci n...n e*], 
where C^(G)r\ e\ n e2 n... ek'-i n (ek') is interpreted to mean all near-min cuts 

containing ei through Ck'-i but not cr- The cuts in this partition, except for the unique cut 
of the last term which has already been found as Co, are enumerated by calling 
ENUMERATE recursively with the argument sets Ef and ET, where Ef denotes the set of 
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Algorithm A1 

INPUT: A directed network G = (V, E), s,t,yv and 8. 

OUTPUT: All minimal s-t cuts C such that w(C) < (1 + e) x w(Co) 

where Co is a minimum weight s-t cut of G. 

begin 

1 Find a min cut Co in G and let Wmin w(Co); 

2 Wmax ^ (1 + £) X Wmin J 

3 <— 0; /* set of edges to be included */ 

4 E~<—0; /* set of edges to be excluded */ 

5 ENUMERATE {G, s, t, w, E\ E', Wmax) ; 

end 

Procedure ENUMERATE {G, s, t, w, E\ E~, Wmax) 

begin 

6 if (3 a cut C' such that E^ cz C' and E~f\ C'—0) 

7 then let C^ Cq and Wnew w( Cq) where Cq is a min 

cut satisfying the conditions for C '; 

8 else return; 

9 if {Wnew > Wmax) then return; 

10 print {Wnew, C^ ); 

11 for (each edge e e { C^ \ £”*' }) 

12 E ^— E o 

13 ENUMERATE {G, s, t, w, , Wmax)', 

14 E~<— E~\ {e}; 

15 E*<^E^u{e}; 

endfor 

16 return; 
end 


Figure 1. Algorithm Al: A generic partitioning algorithm for enumerating all near¬ 
minimum s-t cuts with a weight no greater than 1 + E times the minimum cut weight. 
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edges to be included and ET denotes the set of edges to be excluded. (This essentially 
results in calling ENUMERATE with {(ei)}, {ei, (e 2 )}, {e\, ez, (e 3 )},...,{ci, ez, (ek)} 
for the just-identified cut {ei, ^ 2 , where, as above, edges in parentheses are forced 

out of possible solutions while the others are forced in.) The procedure calls itself 
recursively for every “unforced edge” of the locally minimium cut that is identified 
within each call to ENUMERATE, unless it is determined that no acceptable cuts remain 
and the algorithm must backtrack. 



Figure 2. Sample Graph. Numbers represent edge weights. 


To illustrate how the enumeration works, consider the enumeration tree (which is 
an instance of a rooted tree) given in Figure 3, which corresponds to solving ANMCP-st 
on the graph of Figure 2. The algorithm first finds a minimum cut, {ez, ^ 4 , es) at the root 
node at level 0, then recursively partitions the solution space via {(^ 2 )}, {^ 2 , (^ 4 )} and 
{ez, €4, (cg)}. Once an edge of a cut at some node k has been processed, then it will never 
be processed again at any descendant node of k, because its status as “included” or 
“excluded“ with respect to the current cut has been fixed at node k. The branches with 
{(^ 2 )}, {^ 2 , (^ 4 )} and {ez, €4, (cg)} correspond to searches for a new minimum cut by 
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processing the edges as described. If a search is successful, it leads to a productive node 
where the edges of a new cut, whose status has not been fixed, are processed. Otherwise, 
the search leads directly to a terminal node, along an unproductive branch, where the 
algorithm backtracks. The correctness of this partitioning scheme follows from the fact 
that w(Co) is a monotonic, non-decreasing function of edge weights in G and the 
partitioning is based on a straightforward inclusion-exclusion technique. The actual 
implementation of the algorithm could be difficult, and efficiency poor, because edge 
inclusion may not be so simple to ensure (although exclusion is). 



(Level 0) 


(Level 1) 


(Level 2) 


(Level 3) 


Figure 3. Enumeration Tree for the Graph of Figure 2. The enumeration scheme is 
represented from top to bottom and right to left. Each black node corresponds to a cut 
(edges given in bold next to the node) whose weight is no larger than 1+8 = 1+0.4 times 
the min-cut weight, i.e., Wmax = L(l+0.4)x3j = 4. Numbers in parentheses represent the 
number of the edge excluded from a cut; whereas other numbers represents edges to be 
included. White nodes are terminal nodes where the algorithm backtracks. The 
branches arriving into a terminal node are called “unproductive branches.” 
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We next present a “relaxed” version of Algorithm Al, denoted “Algorithm A2,” 
which may waste time working with non-minimal cuts: It partitions the space of minimal 
and non-minimal cuts but only prints out the minimal cuts. Our final algorithm, 
Algorithm B (Figure 5), closely mimics Algorithm A2 (Figure 4), although Algorithm B 
need not enumerate all near-min, non-minimal cuts. 


Algorithm A2 

INPUT: A directed network G = (V, £), 5, r, w and e 

OUTPUT: All minimal s-t cuts C such that w(C) < (1 + £) x w(Co) 

where Co is a minimum weight s-t cut of G. 

begin 

Same as Algorithm Al; 

end 

Procedure ENUMERATE (G, s, t, w, E*, E~, Wmax) 

begin 

6 if (3 a minimal or non-minimal cut C' s.t. £"^c C' and E~ 0 C'=0) 

I then let and Wnew <-w(Cq) where C'q is a min-weight 

minimal or non-minimal cut satisfying the conditions for C '; 

8 else return; 

9 if (Wnew > Wmax) then return; 

10 if (Cj is minimal) then print (Wnew, ); 

II for (each edge ee { \ }) 

12 E <— E Kj {e}; 

13 ENUMERATE (G, s, t, w, , Wmax)', 

14 E~ <r- E~\ {e}; 

15 E^^E^^{e}; 

endfor 

16 return; 
end 


Figure 4. Algorithm A2: A “relaxed” version of Algorithm Al. 


15 



B. AN IMPLEMENTABLE ALGORITHM 

Algori thm B, given in Figure 5, implements a variant of Algorithm A2, which 

may or may not produce non-minimal cuts. We quasi-exclude an edge e = (m,v) from 
possible cuts at subsequent levels of the enumeration tree by simply setting w(e) = »=. (In 
fact, we use a large number z instead of «>, i.e., z > notational 

convenience we use oo.) It is clear that every near-minimum cut in G should have a finite 
weight, and thus e with w(e) = o® cannot be contained in (G) . This means that quasi¬ 
exclusion implements true exclusion. The graph G with edge e quasi-excluded is denoted 
G-e. 

We quasi-include e = (m,v) by adding two additional edges to G as follows: One 
edge ( 5 ,m) connects the source to m, and the other edge (v,0, connects v to the sink; both 
new edges have infinite weights. The graph G with edge e quasi-included is denoted G + 
e. It is clear that if G + e contains a finite-weight cut, all such cuts must contain e = (m,v) 
because, by construction, u must be on the s side of such a cut, v must be on the t side of 
such a cut, and that implies that e is a forward arc in the cut. Actually, we implement 
quasi-inclusion in Algorithm B by temporarily treating m as an additional source and v as 
an additional sink. 

Figure 6 illustrates the quasi-inclusion and -exclusion of edges. G + -E~ will 

denote G with the set of edges quasi-included and the set of edges E~ quasi-excluded. 
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ALGORITHM B 


INPUT: A directed network G = {V, E), s,t,w and 8 

OUTPUT: All minimal s-t cuts C in G such that w(C) < (1 + 8) x w(Co) 

begin 

1 [Wmin ,Co]<r- MAXFLOW (G , 5, t, w); 

2 Wmax <— (1 + 8) X Wmin ; 

3 0; /* set of edges to be included */ 

4 E~ <r-0\ !* set of edges to be excluded */ 

5 ENUMERATE (G, 5, t, w, E\ E' , Wmax) ; 

end 

Procedure ENUMERATE (G, s, t, w, E'*' , E~ , Wmax) 

begin 

6 w' w; 

7 for (each edge e=(u, v) in E~ ) w'iu, v) 4— oo; 

8 for (each edge e=(u, v) in E'^ ) 

9 add artificial edge (s, u) to G and let w'(s, u) 4— oo; 

10 add artificial edge (v, 0 to G and let w'{v, t) 4— qo; 

endfor 

/* w' is now interpreted to include artificial edges */ 

11 [Wmw ,€,]<- MAXFLOW (G , 5, t, w'); 

12 if {Wnew > Wmax ) then return; 

13 if (is minimal) then print {Wmw, ); 

14 for (each edge e e { \ }) 

15 £"'4- E~yj {e}; 

16 ENUMERATE (G, 5, t, w, E\ E' , Wmax); 

17 E 4— E \{e}; 

18 £"4-£"u{e}; 

19 endfor 

20 return; 
end 


Figure 5. Algorithm B: An approximate implementation of Algorithm A2. 
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C. CORRECTNESS OF THE ALGORITHM 

Algorithm A2 would correctly find all near-minimum minimal cuts if it were 
partitioning the set of all minimal cuts in G along with some, rather than all, of the non- 
minimal cuts. That is essentially what Algorithm B does. Algorithm B begins by finding 
a minimum cut in G and determining Wmax using a max-flow algorithm, all in an 
obviously correct manner. That is, the main routine of Algorithm B correctly implements 
the main routine of Algorithm A2. Then, where Algorithm A2 finds a min cut, minimal 
or non-minimal, that includes edges in £* and excludes edges in E~, Algorithm B solves a 
max-flow problem and finds a min cut, minimal or non-minimal, in G + - E~. 

Algorithm B is clearly finite and will be correct as long as (a) the weight of the locally 
minimum cut is monotonically non-decreasing, and (b) it correctly partitions the space of 
all minimal and possibly some of the non-minimal cuts in G. We first show point (a). 
Lemma 3.1 w(Co) is monotonically non-decreasing for each call to ENUMERATE in 
Algorithm B. 

Proof: When ENUMERATE is called, the finite weight of one edge in G has been 
increased to infinity, and some infinite-weight edges may have been added to G. Clearly, 
this cannot decrease the maximum flow (which equals the minimum cut weight) in G and 
thus w(Co) is monotonically non-decreasing. ■ 

Now, point (b) above will be satisfied if no non-minimal cuts are lost in the calls 
to ENUMERATE and none are repeated. The following two lemmas suffice to prove this. 
Lemma 3.2 Let C be a finite-weight set of edges in G and let E^ and E~ be quasi¬ 
inclusion and quasi-exclusion sets, respectively, produced while running Algorithm B. 
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Suppose that C n iE^ = and C n = 0. Then C is a minimal cut of G only if C is 
also a finite-weight minimal cut of G + - E~. 

Proof: So, we suppose that C is a minimal cut of G. The minimality of this cut is not 
affected by the weight of any edge, so C is also minimal cut of G — E~. C is also a cut of 
(G -E~) + Ef = G -ET + Ef because the only way to stop this, while transforming 
G -ET into G -ET +Ef, would be to add one or more edges from the s side of C to the t 
side of C. But, quasi-inclusion of Ef only adds edges between nodes on the same side of 
the cut. So, C is definitely a cut G - +Ef. Is it a minimal cut? No paths in G - 
have been deleted by adding the edges associated with quasi-inclusion of Ef (to obtain 
G - ET +Ef from G -ET). Thus, all edges of C must be deleted from (G -Er) + E* = 
G — ET + Ef in order to disconnect s from t in that graph; this means that C is a minimal 
cut 'm G - ET + Ef. Finally, C must be a finite-weight minimal cut in G - IT + Ef 
because it contains no infinite-weight edges. ■ 

Lemma 3.3: Let C be a set of edges in G and let Ef and ET be quasi-inclusion and quasi¬ 
exclusion sets, respectively, produced while running Algorithm B. Suppose that 
CnEf ^Ef or CriEr^0. Then, C is a not a finite-weight minimal cut of G + - ET. 

Proof. From a previous argument, we know that quasi-exclusion properly implements 
edge exclusion. Thus, C cannot be a finite-weight minimal cut of G + - F“ if some 

edge of C has been excluded, i.e., ii C n ET ^ 0. Also from a previous argument, we 
know that all finite-weight minimal cuts of G + F*" must contain Ef, i.e., C cannot be a 
finite-weight minimal cut of G + Ef —ET ii C r\Ef ^ Ef. ■ 
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Theorem 3.1 Algorithm B correctly solves ANMCP-st. 

Proof: Lemmas 3.1 through 3.3 show that Algorithm B correctly partitions the set of 
minimal cuts in G (along with, possibly, some non-minimal cuts): The partitioning 
process never loses a minimal cut by Lemma 3.2. Also, the process never repeats a 
minimal cut by Lemma 3.3 and the fact that the algorithm backtracks if the “min-weight 
cut” has infinite weight. Lemma 3.1 shows that the weight of the locally minimum cut in 
Algorithm B is monotonically non-decreasing and, therefore, the algorithm is finite and 
cannot miss enumerating a minimal cut because certain cut weights are “skipped” due to 
backtracking ■ 

D. COMPLEXITY OF THE ALGORITHM 

Algorithm B can be used to solve both AMCP-st and ANMCP-5-r, depending on 

the choice of e. We first give a complexity analysis for all min-cut enumeration, AMCP- 
st. 

1. Complexity Analysis of Minimum-Cut Enumeration 

Consider the enumeration tree given in Figure 3. Every node in that tree is either 

productive, and defines a new cut, or it is a terminal node and is immediately backtracked 

from. Unfortunately, the quasi-inclusion technique might result in the identification of 

one or more non-minimal cuts, i.e., some productive nodes might correspond to non- 

minimal cuts. Fortunately, any non-minimal cut encountered while solving AMCP-st 

cannot be productive and an efficient procedure results. 

We know that the worst-case complexity of solving, from scratch, an initial max 

flow problem on G = (K-E) is 0(f (n,m)), where /(n,m) is a polynomial function of 

n = IVI and m = lEI. At each non-root node of the enumeration tree, the local max flow is 
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obtained by attempting to perform flow augmentations starting with the feasible flow 
from the parent node. (This will be shown to be correct later.) Each flow augmentation 
requires 0(m) work using breadth-first search in a standard fashion. It then turns out that 
the total amount of work performed at each node is 0(m), because (a) if the first search 
does not find a flow-augmenting path, a new min cut has been identified, and (b) if a 
flow-augmenting path is found, the locally maximum flow is at least Wmin + 1 and the 
algorithm can backtrack immediately. (The algorithm must be modified slightly to do 
this.) Now a non-minimal cut must have a weight of at least Wmin + 1, so no such cuts 
correspond to productive nodes. Thus the number of productive nodes is \Cq(G) I. 

Now, each productive node can generate at most n non-productive nodes 
(assuming G has no parallel arcs), so the total number of nodes generated is bounded by 
(n + 1) I Co(G) I. As described above, the amount of work to generate each node except 

the first is 0(m), and the amount of work to generate the first node is 0(f{n,m )), so we 
have the following result. 

Theorem 3.3. Algorithm B finds all minimum cuts (solves AMCP-st) in 
G(/(«,m)-l-/nn ICo(G) I)) time. ■ 

This shows that the time complexity of Algorithm B is polynomial in the number 
of minimum cuts enumerated and the size of graph. Algorithm B is somewhat less 
efficient for solving AMCP-st than are some other algorithms from the literature: The 
algorithm of Picard and Queyranne (1980) solves ANMCP-st in 
G(/(n,m) + (m-i-n)(ICo(G)l)) time and the algorithm of Ball and Provan (1983) solves 

the problem in 0(nnz-i-m I Cg (G) I) time. Nevertheless, our algorithm has several 
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advantages in that (a) it is easy to implement, (b) empirically, it tends to run in 
0(/(«,m) + mlCo(G)l) time and, (c) it extends trivially to near-min cut enumeration, 

i.e., to solving ANMCP-st. 

2. Complexity Analysis of Near-Minimum Cut Enumeration 

By the arguments of the preceding section, the number of productive nodes in 

enumeration tree should be bounded by (« + l)(IC£(G)l + IC'(G)l), where C'(G) 

denotes the set of near-minimum, non-minimal cuts identified as productive nodes of 
Algorithm B’s enumeration tree. The test for non-minimality, which utilizes a breadth- 
first search, takes 0{m) time at each node. The search for a local max flow might require 
multiple flow augmentations and might be as hard as solving for a max flow from 
scratch. Therefore, the work expended at every node is 0{f{n,m) + m) = 0{f(n,m )). If 
we could backtrack whenever a non-minimal cut was identified, then we could say 
IC'(G)I< nlC^(G)l because no node could have more than n descendants and the 

resulting complexity for the whole algorithm would be G(n/(n,m)IC^(G)l). In near¬ 
minimum cut enumeration, however, non-minimal cuts are acceptable if their weight 
does not exceed (1+£)Wmin, and empirical tests show that backtracking when a non- 
minimal cut is encountered can result in the loss of some valid minimal cuts. Thus, 
Algorithm B must ignore non-minimal cuts and continue until it can backtrack based on 
cut weight. This results in a complexity of 0(nf(n,m){\ C^(G) I +1 C'(G) I)), which may 

not be polynomial if I C'(G) lis exponentially larger thanlC^(G) I. Therefore, the worst- 

case complexity of Algorithm B for near-minimum cut enumeration is not well 
determined. We leave this complexity issue as a topic for future research. 
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IV. COMPUTATIONAL RESULTS 


In this chapter, we report computational experiments with Algorithm B to 
demonstrate its correctness and efficiency for solving both AMCP-st and ANMCP-st. 
We test Algorithm B on both weighted and unweighted grid graphs and on various 
problem instances from the literature. 

Algorithm B is written and compiled using the Java 1.2.2 (Sun Microsystems 
1998) programming language. All tests are performed on a personal computer with a 733 
MHz Pentium III processor and 128 MB of RAM, running under the Windows 98 SE 
operating system. 

A. IMPLEMENTATION DETAILS 

1. Efficient Implementation of Algorithm B 

Here, we discuss some improvements of Algorithm B that empirically yield a 
more efficient algorithm for solving both AMCP-st and ANMCP-st. 

The MAXFLOW routine of Algorithm B is implemented to solve an “incremental” 
max flow problem rather than solving for a max flow “from scratch.” This can be done 
because the max flow F* in G is a feasible flow F in G -i- Consider the 

enumeration tree introduced in Chapter III, and let G* and F** be the modified graph and 
its corresponding max flow value, respectively, at node k. Each max flow problem at any 
descendant node of ^ is a relaxation (with respect to edge weights) of the max flow 
problem at node k. Therefore, F/* in Gi, at any immediate child of k, can be reoptimized 
from Ft* in Gt by using an augmenting-path algorithm. (See Ahuja et al., 1993, pp. 180- 
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184 for details.) A complete max flow problem is only solved at root r; the remaining 
max flow problems at any node of the enumeration tree can be obtained by using F* of 
parent nodes. In fact, it takes at most 0{m) work at each node when solving AMCP-st, 
because we need to attempt only a single flow augmentation before we find a new max 
flow or discover that the local max flow is greater than Wmax which is equal to Wmin in 
AMCP-st. Unfortunately, for ANMCP-st, multiple flow augmentations might be needed 
to re-optimize the flow and this might require as much work as solving a complete max 
flow problem. In practice, much less work is required, however. 

The test for the non-minimality of a cut identified in Algorithm B is 

performed only while solving ANMCP-st. The algorithm checks for non-minimal cuts 
right after the identification of each cut in MAXFLOW. It performs a (breadth-first) 
search starting at s trying to reach as many vertices as possible without traversing any cut 
edges and keeps track of the vertices that have been reached. Then, a similar search is 
performed backward from t. A cut is identified as minimal if for every edge 
e = (m,v) e C^, u is reached from 5 and v is reached from t. Otherwise, the cut is non- 

minimal and the algorithm does not print it. 

One other issue in efficient implementation is edge inclusion. In theory, we 
quasi-include an edge (m,v) by adding infinite-weight edges {s,u) and (v,0 to the graph, 
but in practice we simulate this by simply treating m as an additional source and v as an 
additional sink. 

The rest of the implementation is straightforward. We use forward and reverse 
star representation as our data structure (Ahuja et al. 1993, pp. 35-38) and the shortest 
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flow-augmenting-path algorithm of Edmonds and Karp (1972) for solving max flow 
problems. 

2. Problem Generators 

In our literature search, we have not seen any particular problem family designed 
to generate instances for testing algorithms for AMCP-st and ANMCP-st, except for Grid 
Graph Families (GGFs) (Curet 1999, Gibbons 2000). The problem generators, however, 
for MCP, MCP-st (e.g., Levine 1997) and the max flow problems at DIM ACS (The 
Center for Discrete Mathematics and Theoretical Computer, DIMACS 1991) are easy to 
modify for our purposes. Therefore, we use GGF problems and several problem classes 
from Levine (1997) and DIMACS to test Algorithm B. 

We code a GGF Generator (GGFGEN) in Java to generate grid graphs. The 
width W of the graph, and its length L, determine the size of the generated graph. Other 
parameters are e and q, where 8 determines the near-minimum weight criterion and q 
indicates whether the graph is weighted {q = 1) or unweighted (^ = 0). For both weighted 
and unweighted graphs, the edges beginning and ending in s and t, respectively, have 
infinite weights. Every grid vertex u (other than s and t) is connected to each adjacent 
(vertically and horizontally, assuming it exists) vertex v with two directed edges, (m,v) 
and (v,m). Edge weights are 1 in the unweighted case, and pseudo-random, uniformly 
distributed integer weights in range [1,10] otherwise. GGFGEN produces a (connected) 
directed graph with (WL+2) vertices and 2(W+(2WL—W-L)) edges. Figure 7 shows a grid 
graph generated by GGFGEN with inputs W=3,L = 4,q = 0 and any £ > 0. 
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Figure 7. An unweighted, directed grid graph generated by GGFGEN with inputs W =3, 
I = 4, ^ = 0 and any 8 > 0. Edges incident to 5 and t have infinite weights. Other weights 
are all 1. Bi-directional edges between u and v represent two directed edges, (u,v) and 
(v,m). 


We consider different weighted and unweighted instances of GGF with W = L 
(GGF-square) and W« L (GGF-long) for unweighted graphs. Table 1 gives the data for 
the generated problem graphs. 


Problem Name 

W 

L 

E 

q 

GGF-square 

5,10,15,20,25,30, 

40,50,60,70,80,90 

5,10,15,20,25,30, 

40,50,60,70,80,90 

0.0, .05, 
.10,.15 

0,1 

GGF-long 

25 

100,125,150,175, 

200,225,250 

0.0 

0 


Table 1. Problem groups for GGF. A WxL grid of vertices is generated and adjacent 
vertices are connected by bi-directional, vertical and horizontal edges. For q = 0, the 
graph is unweighted (all edge weights are 1), and for q — I, the graph is weighted with 
pseudo-random integer weights in the range [1,10]. 
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We have also chosen two other generators from the literature, implemented in the 
C language and available via Internet for research use. The first is the Double-Cycle 
Generator (DBLCYCLEGEN) (Levine 1997). The single input parameter for 
DBLCYCLEGEN is n, the number of vertices. DBLCYCLEGEN generates two 
interleaved cycles such that the outer cycle has n vertices with the edge weights of 1000 
and 997, and the inner cycle connects every third vertex of the outer cycle with the edges 
of weights 1 or 4. A minimum cut lies in the middle of the graph with a weight of 2000 
and there are many near-min cuts of the weight 2006. 

The second generator is the Acyclic Dense (AD) graph generator from DIMACS 
(1991). AD takes n as its input parameter and generates a fully dense, directed acyclic 
graph with n vertices and m = n (n-l)/2 edges. We replace the pseudo-randomly 
generated edge weights in AD with unit weights to observe the behavior of our algorithm 
on the underlying topologic structure. 

Table 2 gives the generated problem types for DBLCYCLEGEN and AD. 


Problem Name 

n 

£ 

DBLCYC-I 
(Levine 1997) 

500 

0.00,0.10,1.25,1.50, 1.75,2.00 

AD-I 

(DIMACS 1991) 

50 

0.10, 0.20,0.30, 0.40, 

0.50, 0.60,0.70 


Table 2. Problem types for DBLCYCLEGEN and AD. Graphs of type DBLCYC-I are 
generated with DBLCYCLEGEN and have n = 500 vertices and m = 1000 edges on two, 
interleaved, directed cycles. The outer cycle has edge weights of 1000 and 997, and the 
inner one has weights of 1 and 4. The edges between the cycles are weighted so as to hide 
a min-cut with a weight of 2000 in the middle of the graph. Graphs of type AD-I are 
fully dense, directed acyclic graphs with n = 50 vertices and m = 1225 unweighted 
edges. 
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B. THE EXPERIMENTS 

Run times are specified as the total of the time (a) spent in MAXFLOW (MF 
TIME) including the initial max flow at the root node, (b) the time for identifying and 
searching for a new cut (ID TIME), and (c) the time for determining if the new cut is 
minimal (CHECK TIME). All times are given in seconds. Graph input time is short, and 
not reported. 

1. Experiments on Unweighted Graphs 

Table 3 presents mn times of Algorithm B on GGF instances for solving AMCP- 
st. It takes less than 1 second for Algorithm B to identify all minimum cuts in grid 
graphs with sizes up to 402 vertices and 1,560 edges. The number of calls to MAXFLOW 
increases about linearly with n. Because no non-minimal cuts are enumerated, most of 
the running time is spent identifying and searching for a new cut. Also, note that 
CHECK TIME increases as graph size increases. This is due to the 0{m) running time of 
breadth-first search used to check for cut minimality. 
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Problem 

Name 

IVI 

\E\ 

ICol 

ICo(G)l 

Non- 

min. 

Cuts^ 

Calls 

to 

MF^ 

MF 

TIME^ 

(sec.) 


ID 

TIME® 

(sec.) 

TOTAL 

TIME* 

(sec.) 

GGFlOxlO 

102 

380 

10 

9 

0 

92 

0.05 


■EES 

■keei 

GGF20X20 

402 

1560 

20 

19 

0 

382 

0.00 


■GESI 

bees 

GGF30X30 

902 



29 

0 

BBS 

IBESI 

l^EESI 

mm 

^besei 

GGF40X40 

1602 

MeVitl 

40 

39 

0 

Btfgfel 

■El 


^EESi 


GGF50X50 


9900 

50 

49 

0 

2452 

0.11 

0.00 


IKS^i 

GGF60X60 



mi 

59 

0 

lOSEi 

0.21 




GGF70X70 

msmm 

m 

El 

69 



0.48 



■KS 

GGF80X80 

6402 

25440 

80 

79 

0 

6322 

1.15 

HElQi 


K^S 










■■I 


GGF25X100 

2502 

9800 

25 

99 

0 



0.12 


■ESEI 

GGF25X125 

3127 

12250 

25 

124 

0 

3102 

0.06 

0.33 



GGF25X150 

3752 

14700 

25 

■Q3 



0.06 

0.61 


mEESi 

GGF25X175 

4377 

17150 

25 

■■Bl 

0 

1^3 

0.11 

0.90 

4.64 


GGF25X200 

5002 

19600 

25 

199 

0 

4977 

0.11 

1.55 

7.34 


GGF25X225 

5627 


25 



5602 

Kill 

1.88 



GGF25X250 

6252 


25 


0 

6227 


2.13 




1 Non-Min. Cuts : Number of the non-minimal cuts enumerated. 

2 Calls to MF : Number of calls to MAXFLOW . 

3 MF TIME : Total amount of the time spent in MAXFLOW including the initial complete max flow problem. 

4 CHECK TIME : Total time for determining the if each new cut is minimal. 

5 ID TIME : Total time for identifying and searching for a new cut. 

6 TOTAL TIME : Total running time in seconds. 

Table 3. Run times (in seconds) for Algorithm B solving AMCP-st on unweighted 
instances of GGF-square and GGF-long graphs. No non-minimal cuts are enumerated 
because the weight of a non-minimal cut is at least Wmin -t- 1, and the algorithm 
backtracks immediately if a cut with weight greater than Wmin is found. GGFWxL 
denotes a GGF graph with a VFxL grid of vertices. 


Table 4 summarizes the results for ANMCP-st on GGF-square instances with 
£ = 0.05, 0.10, and 0.15. Solution times are expected to increase as £ increases, because 
the number of cuts in any graph might be exponential in the size of the graph. The 
algorithm is quite efficient for modest-size grid graphs with reasonable £ values. 
Compared with Gibbons’ results for near-minimum cut enumeration (Gibbons 2000), our 
results indicate an exponential decrease in calls to MAXFLOW and run times. 
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Problem 

Name 

IVI 

\E\ 

ICol 

ICJ 

ICe(G)l 

Non- 

min. 

Cuts 

Calls to 
MF 

MF 

TIME 

(sec.) 

CHECK 

TIME 

(sec.) 

ID 

TIME 

(sec.) 

TOTAL 

TIME 

(sec.) 

e = 0.05 



GGF5x5 

27 

90 

5 

5 

4 

0 

22 

0.05 

0.00 

0.05 

O.IO 

GGFlOxlO 

102 

380 

10 

10 

9 

0 

92 

0.00 

0.00 

0.06 

0.06 

GGF15X15 

227 

870 

15 

15 

14 

0 

212 

0.00 

0.06 

0.00 

0.06 

GGF20X20 

402 

1560 

20 

21 

703 

0 

7906 

0.77 

0.00 

0.71 

1.48 

GGF25X25 

627 

2450 

25 

26 

1128 

0 

15506 

1.62 

0.15 

2.14 

3.91 

GGF30X30 

902 

3540 

30 

31 

1653 

0 

26856 

3.92 

0.49 

7.55 

11.96 

e = 0.10 

GGF5x5 

27 

90 

5 

5 

22 

0 

22 

0.05 

0.00 

0.05 

0.10 

GGFlOxlO 

102 

380 

10 

11 

956 

0 

956 

0.00 

0.06 

0.05 

0.11 

GGF15X15 

227 

870 

15 

16 

3306 

0 

3306 

0.23 

0.00 

0.27 

0.50 

GGF20X20 

402 

1560 

20 

22 

113090 

0 

113090 

9.11 

1.35 

9.63 

20.09 

GGF25X25 

627 

2450 

25 

27 

274550 

0 

274550 

30.98 

3.83 

39.62 

74.43 

GGF30X30 

902 

3540 

30 

33 

8911698 

378 

8911698 

1553.89 

182.78 

2798.42 

4535.09 

e = 0.15 


GGF5x5 

27 

90 

5 

5 

22 

0 

22 

0.05 

0.00 

0.05 

0.10 

GGFlOxlO 

102 

380 

10 

11 

956 

0 

956 

0.05 

0.05 

0.06 

0.16 

GGF15X15 

227 

870 

15 

17 

35905 

0 

35905 

1.69 

0.42 

2.77 

4.88 

GGF20X20 

402 

1560 

20 

23 

1202033 

153 

1202033 

102.92 

13.70 

98.70 

215.32 

GGF25X25 

627 

2450 

25 

28 

3621978 

253 

3621978 

420.08 

56.83 

496.38 

973.29 

GGF30X30 

902 

3540 

30 

34 

13465371 

21843 

113463496 

18537.49 

2609.65 

25248.65 

46395.79 


Table 4. Run times in seconds for ANMCP-st solutions of Algorithm B on unweighted 
GGF-square instances with e = 0.05, 0.10, 0.15. As expected, the solution times increase 
as £ increases. For example, 8,911,698 near-min cuts together with 378 non-mioimal cuts 
are enumerated in 4,535.09 seconds in GGF30x30 for E = 0.10. For £ = 0.15 , on the 
same graph, 13,465,371 near-min cuts and 21,843 non-minimal cuts are enumerated in 
46,395.79 seconds. 


Table 5 presents results for unweighted AD graphs. Interestingly, Algorithm B 
does not enumerate any non-minimal cuts even though £ is relatively large in some 
instances. Given our concern about enumerating non-minimal cuts, this invites further 
investigation. 
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e 

ICol 

\CJ 

IQCG)! 

Non- 

min. 

Cuts 

Calls to 
MF 

MF 

TIME 

(sec.) 

CHECK 

TIME 

(sec.) 

ID 

TIME 

(sec.) 

TOTAL 

TIME 

(sec.) 

0.0 

49 

49 

49 

0 

1275 


0.00 

0.11 

0.32 

0.10 

49 

53 

544 

0 

13650 

2.28 

0.00 

0.85 

3.13 



58 

4063 

0 

101625 

16.56 

0.22 

8.92 


MiKM 

49 

63 

19798 

0 

495000 

94.86 

0.64 

39.40 


0.40 

49 

68 

75893 

0 

1897375 

388.38 

3.96 

150.61 

542.95 

0.50 

49 

73 

249270 

0 

6231800 

1356.31 

15.90 

489.37 

1861.58 

0.60 

49 

78 

730603 

0 

18265125 

4077.99 

42.87 

1423.98 

5544.84 

0.70 

49 

83 

1962849 

0 

49071275 

11188.42 

121.20 

3816.87 

1 - 

15126.49 


Table 5. Computational results on unweighted, acyclic dense (AD) graphs with various 
threshold levels £. This is a fully dense graph with 50 vertices and 1225 edges. Although 
the graph is fully dense, the number of non-minimal cuts enumerated for all values of e is 
zero. 


2. Experiments on Weighted Graphs 

Here we use the GGF-square problems with edge weights pseudo-randomly 
generated integers in the range [1,10]. Results for minimum and near-minimum cut 
enumeration are summarized in Tables 6 and Table 7, respectively. 


Problem Name 

IVl 

l£l 


iCo(G)l 




CHECK 

TIME 

(sec.) 


TOTAL 

TIME 

(sec.) 

GGF5x5w 

27 

90 

17 

1 

0 

7 

0.06 

0.00 

0.06 


GGFlOxlOw 

102 



2 

0 

92 

0.00 

0.00 

0.05 


GGF15xl5w 

227 


wm 

1 

0 

19 

0.00 

0.00 

0.05 


GGF20x20w 

402 

1560 

69 

1 


32 

0.06 


0.06 


GGF25x25w 

627 

2450 


1 

HH 

33 



0.05 


GGF30x30w 

902 

3540 

108 

6 

0 

217 

0.05 

0.06 

0.16 

0.27 


Table 6. Computational results for minimum-cut enumeration (AMCP-st) on weighted, 
GGF-square problems. As in the unweighted case, no non-minimal cut are encountered. 
All minimum cuts are identified in less than one second for these instances. 
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Problem 

Name 

\V\ 

\E\ 

^MIN 

^MAX 

ICe(G)l 

Non- 

min. 

Cuts 

Calls to 
MF 

MF 

TIME 

(sec.) 

CHECK 

TIME 

(sec.) 

ID 

TIME 

(sec.) 

TOTAL 

TIME 

(sec.) 

£ = 0.05 




GGF5x5w 

27 

90 

17 

17 

1 

0 

7 

0.05 

0.00 

0.05 

0.10 

GGFlOxlOw 

102 

380 

42 

44 

5 

0 

46 

0.00 

0.00 

0.05 

0.05 

GGF15xl5w 

227 

870 

45 

47 

3 

0 

43 

0.00 

0.05 

0.00 

0.05 

GGF20x20w 

402 

1560 

69 

72 

20 

0 

308 

0.11 

0.00 

0.06 

0.17 

GGF25x25w 

627 

2450 

87 

91 

33 

0 

645 

0.16 

0.00 

0.11 

0.27 

GGF30X30W 

902 

3540 

108 

113 

416 

0 

7221 

1.59 

0.17 

2.04 

3.80" 

e = 0.10 



GGF5x5w 

27 

90 

17 

18 

1 

0 

7 

0.06 

0.00 

0.06 

0.12 

GGFlOxlOw 

102 

380 

42 

46 

11 

0 

78 

0.05 

0.00 

0.05 

0.10 

GGF15xl5w 

227 

870 

45 

49 

6 

0 

75 

0.00 

0.00 

0.06 

0.06 

GGF20x20w 

402 

1560 

69 

75 

116 

0 

1508 

0.22 

0.00 

0.16 

0.38 

GGF25x25w 

627 

2450 

87 

95 

309 

0 

4263 

0.51 

0.11 

0.69 

1.31 

GGF30x30w 

902 

3540 

108 

118 

8025 

0 

104836 

25.60 

1.85 

31.66 

59.11 

e = 0.15 



GGF5x5w 

27 

90 

17 

19 

1 

0 

7 

0.05 

0.00 

0.05 

0.10 

GGFlOxlOw 

102 

380 

42 

48 

22 

0 

161 

0.00 

0.06 

0.00 

0.06 

GGF15xl5w 

227 

870 

45 

51 

22 

0 

223 

0.00 

0.06 

0.00 

0.06 

GGF20x20w 

402 

1560 

69 

79 

1031 

0 

10744 

1.54 

0.06 

0.88 

2.48 

GGF25x25w 

627 

2450 

87 

100 

3345 

22 

35917 

5.11 

0.55 

5.49 

11.15 

GGF30x30vt; 

902 

3540 

108 

124 

144511 

0 

1503784 

348.56 

32.41 

466.98 

847.95 


Table 7. Computational results for near-minimum cut enumeration (ANMCP-st) on 
weighted GGF-square problems. Non-minimal cuts are only encountered GGF25x25w. 


Finally, we test Algorithm B on the DBLCYC-I problems with e values ranging 
from 0.0 to 2.0. These tests are the only problem types where Algorithm B generates 
large numbers of non-minimal cuts compared to the number of near-minimum cuts. For 
e > 1.25, the ratio of the number of non-minimal cuts to the number of minimal cuts 
increases dramatically; see Table 8. 
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e 


^MAX 

ICe(G)i 

Non-min. 

Cuts 

Calls to 
MF 

MF 

TIME 

(sec.) 

CHECK 

TIME 

(sec.) 

ID 

TIME 

(sec.) 

TOTAL 

TIME 

(sec.) 


1000 

1000 

2 

HIHQ] 

10 

HEESI 




HQQHI 

1000 

1100 

499 


1976 



0.28 





511 

8 


■ebi 


■01^ 


1.25 


MSSSl 








1.50 

liigg 



237411 

957178 

KkUl 




1.75 

1000 

2750 

2479 

237411 

957178 

92.76 

23.16 

91.48 

201 AO 

2.00 

1000 

3000 

2509 

238041 

959683 

89.67 

21.01 

82,72 

193.40 


Table 8. Computational results for near-minimum cut enumeration (ANMCP-st) on 
weighted DBLCYC-I problems when n = 500 (and m = 1000). The number of the non- 
minimal cuts increases dramatically with increases in e. For example, when 8 = 2.00, i.e., 
Wmax = (l+2.00)x 1000 = 3000, the number of non-minimal cuts is 238,041 whereas the 
number of near-min cuts is only 2509. 


Our results show that Algorithm B performs quite well on various types of graphs, 
but that non-minimal cuts can slow computations when the threshold parameter e 
becomes large. For dense acyclic graphs, Algorithm B does not enumerate any non- 
minimal cuts, but for the “double-cycle graphs” DBLCYC-I, the number of non-minimal 
cuts can outnumber the minimal cuts by a huge margin, at least when 8 becomes large. 
This behavior of the algorithm suggests that it does have polynomial complexity for 
certain graph topologies, but not others. We believe this issue merits further 
investigation. 

For large graphs where the time for flow augmentations dominates the overall 
running time, use of a state-of-the-art max flow algorithm such as Goldberg and Rao 
(1998) might reduce MAXFLOW time. 
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V. CONCLUSIONS AND RECOMMEDATIONS 


In this thesis, we have developed an algorithm for enumerating all near-minimum- 
weight s-t cuts in a directed graph G = (V,E). The enumeration algorithm finds a 
minimum-weight cut in the input graph via a maximum flow algorithm and then 
recursively partitions the space of possible solutions to the “near-min threshold level” 
Wmax = (l+e)xWW/Af, where Wmin is the weight of that minimum weight cut. Given a 
minimum cut, the set of acceptable cuts is recursively partitioned by forcing inclusion 
and exclusion of edges from subsequent cuts. An edge (m,v) is quasi-excluded by simply 
setting its weight to the infinity and quasi-included by implicitly introducing two infinite- 
weight edges in G, one extending from ^ to m and the other from v to t. The algorithm 
solves a max flow min-cut problem for each modified graph that is obtained in the 
enumeration tree. 

We have implemented our algorithm using the following enhancements to 
improve the efficiency: (1) The algorithm solves a complete max flow problem at the 
beginning (at the root node of enumeration tree) but solves only an “incremental” max 
flow problem at the all other nodes (the max flow at a parent node is feasible for all 
descendants and can easily be re-optimized), and (2) quasi-inclusion of an edge (m,v) is 
simulated by treating u as an additional source and v as an additional sink. 

Unfortunately, our quasi-inclusion technique sometimes leads to the enumeration 
of non-minimal cuts together with the minimal cuts. Non-minimal cuts are easily 
identified (and ignored), but they can increase the computational workload and stop us 


37 



from deriving a polynomial-time bound for the worst-case complexity of the algorithm. 
We do obtain, however, a polynomial bound of 0{f{n,m) + nm\Cg{G)\) for minimum- 

cut enumeration (which is the same as near-minimum cut enumeration but with e = 0). 

Computational results for e > 0 show that the algorithm is empirically efficient for 
modest-sized problems with modest values for £. For example, in an unweighted grid 
graph with 402 vertices and 1,560 edges, (a) all 19 minimum cuts can be enumerated in 
0.05 seconds (on a 733 MHz Pentium III personal computer), (b) it takes 0.77 seconds to 
enumerate all 703 near-minimum cuts when e = 0.05, and (c) it takes only 20.09 seconds 
to enumerate all 113,090 near-minimum cuts when e = 0.10. However, in another 
unweighted graph with 102 vertices and 380 edges, and e = 0.50, the algorithm 
enumerates 134,705 near-min cuts together with 4,474 non-minimal cuts in 601.2 
seconds. For e = 0.90, on the same graph, 7,811,043 near-min cuts and 1,941,792 non- 
minimal cuts are obtained in 2,303.5 seconds. The running times and the number of non- 
minimal cuts increase significantly for relatively large values of e. 

Although Algorithm B appears to be quite efficient in practice, more work might 
improve the quasi-inclusion technique or develop another technique for edge inclusion. 
If “true edge inclusion” (as opposed to quasi-inclusion) can be efficiently implemented, 
this should yield a provably polynomial-time algorithm for near-minimum cut 
enumeration. 

If the current quasi-inclusion technique is retained, another approach might be 
used to avoid enumerating non-minimal cuts. In particular, edges that cannot occur in 
any minimal cut given those that are already included can be identified and marked as 



“forbidden for inclusion.” This means that they can be given an infinite weight and 
excluded because they are not allowed to appear in any minimal cut given the current 
inclusions. An edge (u,v) can be forbidden from inclusion if every (s-u) or (v-f) path 
contains at least one such included edge. 

Another approach to avoid non-minimal cuts might be to take the advantage of 
max flow min-cut duality and corresponding linear programming (LP) techniques. This 
might be fruitful because our algorithm enumerates some of the near-minimum-cost 
bases in the dual of the maximum-flow LP. Although a minimal cut might correspond to 
more than one basis, a non-minimal cut certainly should be indicated by the absence of a 
basis. Use of linear algebra (linear independence or dependence of basic variables) to 
check the corresponding bases during the enumeration might lead to a new technique to 
eliminate non-minimal cuts. 

Finally, it would be interesting to see if the enumeration algorithm, as it exists 
now, will enumerate only minimal cuts for certain types of graph topologies, e.g., 
undirected s-t planar graphs. Using the dual of a planar graph and shortest-path 
techniques, it is possible to enumerate near-min cuts in an undirected s-t planar graph in 
polynomial time per cut. Thus, it is natural to wonder if our algorithm can also 
enumerate such cuts efficiently. 
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